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ABSTRACT 

We analyze 3'.'5 resolution, high sensitivity radio continuum images of the nearby spiral galaxy NGC 
4258 at 6 and 20 cm derived from multiple observations used to monitor the radio supernova SN 1981K 
(Van Dyk et al. 1992, [ApJ, 396, 1995]). Seven bright H II region and five supernova remnant candidates 
are identified. Extinctions to the H II regions are estimated for the first time from a comparison of radio 
flux densities to new optical fluxes derived from Ha observations by Dutil & Roy (1999, [ApJ, 516, 62]). 
The bright end of the H II region luminosity function is established at each wavelength. The luminosity 
functions are best fit by power laws consistent with the shape of previously published radio and optical 
luminosity functions for a number of galaxies. The supernova remnants are all about 2-3 times the radio 
luminosity of Cas A. In addition, the galaxy's nucleus is shown to have been variable over the SN 1981K 
monitoring period. The spectral index (a) distribution of the anomalous radio arms is investigated and 
found to be relatively uniform at a — —0.65 ± 0.10. 

Subject headings: galaxies: individual (NGC 4258, M106) — radio continuum: galaxies — H II regions 
— supernova remnants — dust, extinction 



1. INTRODUCTION 

Radio continuum emission from normal galaxies is a 
combination of both thermal and nonthermal discrete 
emission from individual sources and diffuse emission from 
the galactic disk. The individual thermal emitters are as- 
sumed to be free-free emitting H II regions, which affect 
the interstellar medium (ISM) through the ionizing radia- 
tion of young stars, and the nonthermal emitters are gener- 
ally thought primarily to be supernova remnants (SNRs), 
which also profoundly affect the ISM by injecting enor- 
mous amounts of kinetic energy and enriched material. 
SNRs are the leading candidates for the production of cos- 
mic rays and contributors to the relativistic ISM. Thus, 
studies of H II regions and SNRs provide important infor- 
mation on the beginning and endpoints of stellar evolution 
and reveal much about the evolution of the ISM and, ul- 
timately, the evolution of galaxies. 

Analyses of extragalactic SNRs and H II regions at ra- 
dio wavelengths are not hampered by the distance ambi- 
guities of Galactic samples and extinction by dust and gas. 
A relatively well-defined samples of SNR and H II region 
distributions in external galaxies can be used to yield an 
unbiased statistical picture of group properties, which can 
yield important information on star formation and the role 
of supernova (SN) and SNR shock waves in determining 
the energetic and dynamical evolution of the ISM, and, 



perhaps, even in triggering new epochs of star formation. 
Radio observations of H II regions are nearly a direct mea- 
sure of the ionizing luminosity of the embedded stars. A 
comparison with optical line luminosities can result in esti- 
mates of the line-of-sight extinction, providing an alterna- 
tive to hydrogen line ratio methods alone, and can make 
it possible to trace the distribution of dust and the na- 
ture of extinction gradients in other galaxies. Extinction 
measurements can also help determine the amount of dust 
internal to individual H II regions, which are essential for 
accurate estimates of atomic and molecular abundances 
and chemical gradients within galaxies. Comparisons be- 
tween thermal radio and far-infrared properties can lead 
to a better understanding of the young stars that power 
H II regions and their role in the galactic-scale, far-infrared 
emission properties. 

Due to limited resolution and sensitivity, large samples 
of discrete radio sources in external galaxies have only been 
observed in the Local Group (e.g., in M31, Braun & Wal- 
terbos 1993; in M33, Viallefond & Goss 1986, Duric et 
al. 1993, Duric et al. 1995, and Gordon et al. 1999; in the 
LMC and SMC, Mills 1983 and Mathewson et al. 1983), 
in intense starburst galaxies (e.g., in M82, Kronberg, Bier- 
man, & Schwab 1985; in NGC 253, Ulvestad & Antonucci 
1997), in M51 (van der Hulst et al. 1988), in M81 (Kauf- 
man et al. 1987), and in NGC 6946 (Lacey, Duric, & Goss 
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1997, Hyman et al. 2000). As part of a larger study of the 
continuum emission from nearby galaxies, we consider in 
this paper the case of NGC 4258. 

NGC 4258 (M106) is a large (18.'6 x 7.'2; NASA/IPAC 
Extragalactic Database, NED), nearby (7.3 Mpc; Her- 
rnstein et al. 1997) spiral galaxy with high inclination 
(i = 71°, Sandage & Tammann 1981), at high galactic 
latitude (b = 68?8), with correspondingly low foreground 
extinction (As = 0.00; NED). The galaxy is also notable 
for its LINER/Seyfert nuclear properties, for the detection 
of a very bright radio continuum nuclear source (Turner & 
Ho 1994), for the discovery of H 2 masers and for the pos- 
sibility of a supermassive black hole at the center (Miyoshi 
et al. 1995; Greenhill et al. 1995). Ha imaging by Courtes 
et al. (1993) and Dutil & Roy (1999) show a number of 
very bright H II regions along the inner spiral arms of the 
galaxy. 

NGC 4258 has been previously studied at numerous 
radio wavelengths, van der Kruit, Oort, & Mathewson 
(1972) observed the galaxy at 20 cm with the Westerbork 
Synthesis Radio Telescope (WSRT) at 24" x 32" resolu- 
tion, revealing very bright and extensive "anomalous" ra- 
dio arms in the east-west direction, de Bruyn (1977), with 
the WSRT at 6 and 49 cm with resolution 7" x 9" and 
56" x 76", respectively, found a nonthermal spectral in- 
dex for the anomalous arms, van Albada & van der Hulst 
(1982), with the WSRT (6"5 beam) and the Very Large 
Array (VLAQ; 11"25 beam), both at 20 cm, found that the 
anomalous arms are multiply branched. Krause, Beck, & 
Klein (1984), with the Effelsberg 100-m telescope at 6.3 
and 2.8 cm (~1' resolution), detected no significant spec- 
tral index variations along the anomalous arms. Hummel, 
Krause, & Lesch (1989) conclude, based on their analysis 
of the polarized emission using the VLA at 6 and 20 cm 
(14" resolution) and Effelsberg at 1.2 cm (40" resolution), 
that the anomalous arms lie in the plane of NGC 4258. 

Although the anomalous arms coincide with faint Ha 
emission (Courtes & Cruvellier 1961; Courtes et al. 1993; 
Dutil & Roy 1999), they were originally thought to be out 
of the plane of the optical galaxy (Courtes & Cruvellier 
1961; van der Kruit et al. 1972), due to the absence of 
H II regions along them and other considerations. The ra- 
dio arms have been modeled as arising from the expulsion 
of gas from the nucleus either in a single event (van der 
Kruit et al. 1972) or by continuous ejection in a jet (Mar- 
tin et al. 1989; Ford et al. 1986; Cecil, Wilson, & Tully 
1992; Pietsch et al. 1994). However, the CO(1-0) obser- 
vations by Krause et al. (1990) revealed a large quantity 
of molecular gas along the radio arms which would have 
required an unreasonably large energy output in the sin- 
gle event model. Recently, Cecil et al. (2000), with VLA 
20 cm observations at 1"3 resolution and high- resolution 
optical data, conclude that the arms arise from changing 
jet activity at the nucleus. Cox & Downes (1996), from 
CO(2-l) observations of the inner parts of the anomalous 
arms, offer an alternative explanation, concluding that the 
arms arise from gas flows and shocks due to a galactic bar 
potential within the disk of NGC 4258. 

The previous radio studies of NGC 4258 focused mostly 
on the analysis of the anomalous arms and not on the dis- 



crete sources in the galaxy, other than the nucleus. We 
present here our high resolution (3"5) and high sensitiv- 
ity (~25 /xJy beam -1 ) observations, which result in a de- 
tailed analysis of the H II regions and SNRs. Although our 
study does not concentrate on the nature of the anomalous 
arms, our radio maps show significantly smaller structures 
in the anomalous arms and smaller-scale spectral index 
variations within them. Optical-to-radio comparisons are 
used to estimate the visual extinctions toward the H II 
regions. Also, the bright end of the H II region luminosity 
function is produced, and the star formation properties of 
the brightest H II regions are also presented. The SNRs are 
shown to be very luminous in the radio, and the nucleus 
appears to exhibit variability over more than a decade. 

2. RADIO OBSERVATIONS 

The Type II radio SN 1981K in NGC 4258 was moni- 
tored with radio continuum VLA observations at 6 and 20 
cm from 1985 March through 1990 May in all four array 
configurations. Details of the data reduction for these ob- 
servations and the results of the analysis for SN 1981K are 
described in Weiler et al. (1986) and Van Dyk et al. (1992). 
We have used the AIPS task DBCON to combine the ob- 
servations at each frequency from 1985 March 30 through 
1990 May 29, and from additional SN 1981K monitoring 
observations made on 1991 July 29, 1991 October 31, and 
1994 May 03 in A and B configurations. The integrated 
times-on-source for the concatenated uv databases are 16.6 
hr at 6 cm and 10.6 hr at 20 cm. The uv data were 
Fourier-transformed and beam-deconvolved using two dif- 
ferent procedures in AIPS. First, the task IMAGR was 
used. Images with equal resolutions of 3" 4 x 3"3 at both 
wavelengths were generated by appropriately weighting 
the data with UVTAPER. These images were used to ana- 
lyze most of the discrete sources in and background sources 
surrounding NGC 4258. However, the anomalous radio 
arms have a very speckled appearance in these images, 
due to cleaning of the very extended emission with a nar- 
row beam. For that reason, the task VTESS, running the 
Maximum Entropy Method (MEM), was used. While the 
CLEAN algorithm in IMAGR is more reliable for imaging 
discrete sources, MEM is more suitable for imaging ex- 
tended emission. Zero-spacing fluxes of 0.40 Jy (Krause et 
al. 1984) and 0.82 Jy (Hummel et al. 1989), corresponding 
to the integrated flux densities at 6 and 20 cm, respec- 
tively, were input to VTESS to correct for missing short 
spacings in the uv coverage. The MEM images of NGC 
4258 are shown in Figures la and b. They have nearly 
the same resolution (3"5 x 3"3) and sensitivity (20 /LtJy 
beam -1 and 30 /xJy beam -1 at 6 and 20 cm, respectively) 
as the images made with IMAGR (not shown) . The labels 
in Figure la refer to discrete sources discussed in §4.1.1 
and §4.1.2 below. 

Corrections for primary beam attenuation have been ap- 
plied to both the CLEAN and MEM images, since the 
anomalous arms extend ~8' east-west and the prominent 
inner north-south spiral arms extend ~6', well within the 
30' FWHM primary beam at 20 cm, but less so, within 
the ~9' (FWHM) primary beam at 6 cm. Since the field 
center of the observations is at the location of SN 1981K 
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(a = 12 h 18 m 59.42 s , 5 = +47°19'31'/0; J2000), the 50 MHz 
bandwidth at both wavelengths did not result in significant 
bandwidth smearing for sources along the northern arm. 
However, smearing is significant at 20 cm for sources in the 
southern arm, where the reduction in the peak flux den- 
sity of a point source is as much as a factor of two (Bridle 
& Schwab 1988). While the integrated flux density of a 
bandwidth-smeared bright source, such as the prominent 
source at the end of the southern spiral arm, should not 
be affected, fainter sources could be smeared into the sur- 
rounding arm emission and noise, and hence, be difficult to 
detect. The detection of faint radio sources in the southern 
arm is discussed further in the next section. 

3. SOURCE DETECTION AND MEASUREMENT 

Because of the complexity of the radio emission from 
NGC 4258, it was often difficult to resolve discrete 
sources from background disk and anomalous arm emis- 
sion. Therefore, only the brightest sources (nos. 1-9; Fig- 
ure la), a few fainter sources in the southern arm (nos. 10- 
12; Figure 3), one background source at the tip of the 
southeastern radio arm, and nine background sources in 
located outside the field of Figures la and b were chosen 
for analysis. The CLEAN maps were used to measure the 
properties of sources 1-9, the nucleus, and the background 
sources. The MEM images were used to measure sources 
10-12 and to generate a spectral index map of the anoma- 
lous radio arms (see §4.4.2). In addition to our radio data, 
Y. Dutil and J.-R. Roy kindly provided to us narrow-band, 
continuum-subtracted Ha imaging data for NGC 4258 (see 
Dutil & Roy 1999 for details). Briefly, the Ha images were 
acquired with the 1.6 m telescope at Mont Megantic Ob- 
servatory, the two satellite [N II] lines were excluded using 
a 10 A-wide filter, the seeing was ~2-3", and the optical 
flux calibration uncertainty is ~20%. We registered the 
radio and Ha images, and show the Ha map overlaid on 
the 20 cm MEM greyscale image in Figure 2. 

The positions of most of the brightest radio sources have 
optical counterparts with peak pixel positions that agree to 
within l'/5. The brightest portions of the anomalous arms 
also coincide with diffuse optical emission, although sig- 
nificantly more structure is apparent in the radio images 
(see §4.4.1). In contrast to the radio images in Figures 
la and b, the optical-radio comparison in Figure 2 shows 
significant Ha emission and many discrete sources in the 
southern arm, besides the very bright source 9. Because of 
the confusion of faint sources with background emission, 
we have used the MEM images to identify sources in the 
southern arm. An MEM-derived contour map of the 20 
cm radio emission from the southern arm and surrounding 
regions is shown in Figure 3. The radio map shows faint 
emission coincident with several of the discrete Ha sources 
in the southern arm. 

The flux densities of the radio sources were measured us- 
ing the same technique as that used for NGC 6946 (Lacey 
et al. 1997; Hyman et al. 2000). First, the total flux was 
measured within an aperture including the source and an 
annulus outside the aperture measuring the background. 
The average background was computed and subtracted 
from the flux within the source aperture to obtain the 
source flux density. Note that source 4 appears to have 
two components in the 20 cm image (Figure la). How- 



ever, since the corresponding source on the 6 cm and Ha 
images does not show any structure and covers the same 
area as the 20 cm source, we have simply measured the 
flux using one inclusive aperture. Fluxes were measured 
for each of the optical counterparts to radio sources using 
the same method. Since the resolution of the radio maps is 
similar to the seeing of the optical observations, identical 
source apertures were used. 

4. DISCUSSION 

4.1. Source Criteria 

The thermal radio emission from H II regions is identi- 
fiable by its flat spectral index, a (where S oc z^ +a ), with 
a nominal value a = —0.1. For nonthermal emission from, 
e.g., SNe, SNRs, and optically thin active galactic nuclei, 
a is steeper, in the approximate range —0.4 to —1.0. Un- 
fortunately, a nonthermal plerionic SNR, which also has a 
flat spectrum (-0.3 < a < -0.1; Weiler & Panagia 1978), 
might be confused with an H II region by this criterion. 
However, extrapolating from the small number of pleri- 
ons compared to the shell-type remnants (Srinivasan et 
al. 1984), it is likely that we would correctly identify the 
SNRs in NGC 4258 down to the surface brightness limit, 
with only one or two plerions possibly misidentified. 

4.1.1. Discrete Sources in NGC 4258 

In Table 1 we present the positions, flux densities, and 
spectral indices for 12 discrete radio sources, including the 
galactic nucleus. The position listed for each source corre- 
sponds to the centroid at 20 cm. SN 1981K is not included 
in the table. Although the uncertainties in the spectral 
indices are appreciable, we employ the criteria above to 
discriminate between thermal and nonthermal sources. Of 
the 13 sources, sources 1-6 and 9 are considered thermal, 
with a between +0.05 and —0.25, and are, therefore, clas- 
sified as H II region candidates. Two bright sources (7 
and 8) and the faint southern arm sources (10, 11, and 
12) have nonthermal a between —0.37 and —0.82, and are 
classified as SNR candidates. 

We have identified and measured optical counterparts 
on the Ha image to the 13 radio sources in Table 1, ex- 
cept for two SNR candidates (sources 7 and 8) and source 
6, which is blended with source 5 on the Ha map (see 
below). The optical fluxes, as well as estimates of visual 
extinction toward the H II regions, are listed in the table 
(see also §4.1.5). The nucleus, source 2, and the nonther- 
mal sources 11 and 12 are coincident with, or are near, 
four X-ray sources (X18, X16, X19, and X14, respectively; 
Vogler & Pietsch 1999). The X-ray counterparts to sources 
11 and 12 are indicated by boxes in Figure 3. 

Note from Figure la that one optical source covers an 
area that includes both radio sources 5 and 6; therefore, a 
separate aperture including both sources was used to mea- 
sure the optical flux and the flux density of the combined 
pair of radio sources (source "5+6" in Table 1). Also, 
while an optical source is located 4'.'2 from the nonther- 
mal source 7 (i.e., somewhat more than one synthesized 
half-power beam width), we have not included its flux in 
the table, since typically the radio and optical source po- 
sitional uncertainties are less than ~1". 
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4.1.2. Background Sources 

Based on estimates of the background source counts 
with flux densities above 100 ^Jy at 21 cm (Oort 1987), 
we expect ~4 background sources to occur in the field. It 
is therefore possible that extragalactic radio sources with 
nonthermal spectra could be confused with our SNR can- 
didates, or that a flat spectrum galaxy nucleus could be 
confused with our H II region candidates. However, since 
all thermal sources and all three southern arm SNR can- 
didates also have Ha counterparts, these sources are less 
likely to be background objects. Although the nonthermal 
source 7 does not have a clear optical counterpart, its lo- 
cation in the northern arm near the central region of the 
galaxy is supportive of an identification as a SNR. The 
nonthermal source 8 could possibly be a background ob- 
ject, since it is has no optical counterpart and is located 
in the interarm region. 

One notable, very bright, nonthermal radio source 
within the field of Figures la and b (source 15), located 
east of the nucleus at the tip of the southeastern radio 
arm, was rejected as a SNR candidate because of its pos- 
sible "head-tail" structure. The 20 cm position of the 
source's bright (bandwidth smeared) "head" component 
is a = 12H9 m 18!92, S = +47°19'07'.'6 (J2000). The bifur- 
cated radio "tail" extends northward. At 6 cm, only the 
northeastern part of the tail is clearly visible above the 
noise. The head has a 20 cm flux density of 2.09±0.23 mJy 
and a nonthermal spectral index a = —0.69 ± 0.15. Hum- 
mel et al. (1989) and van Albada & van dcr Hulst (1982) 
also detected this source in their lower-resolution maps 
and suggested that, given its extreme brightness, nonther- 
mal spectral index, and possible complex structure, it is 
a background radio galaxy rather than a source associ- 
ated with NGC 4258. From its morphology evident in our 
higher-resolution maps, we agree with that conclusion. 

In Table 2 we list source 15, described above, two clus- 
ters of sources (16 and 17), and two additional sources 
(13 and 14), all lying within the 20% attenuation level 
of the primary beam at 6 cm, but, with the exception 
of source 15, outside the field shown in Figures la and 
b. These sources have no optical counterparts. Nine of 
the ten sources listed in the table have nonthermal spec- 
tral indices ranging from a = —0.7 to —1.1. The cluster 
of sources 16a-d corresponds to source 6 in Hummel et 
al. (1989), which they list as a "triple," and the cluster 
of sources 17a-c corresponds to their source 2, listed as a 
"double." Hummel et al. (1989) and van Albada & van der 
Hulst (1982) suggest that the clusters constitute two radio 
galaxies. Indeed, the source cluster 16a~d may consist of a 
flat spectrum (a = — 0.15±0.11) nuclear source (16b) with 
nonthermal radio lobes to the north (16a) and south (16c, 
d). Similarly, the sources 17a and 17c may be the lobes 
of a radio galaxy, with the nucleus located at the position 
of source 17b. While additional sources are visible within 
the larger field-of-view at 20 cm, they are extremely band- 
width smeared; hence, no measurements were attempted 
for these sources, especially since no information is avail- 
able for them at 6 cm. 

4.1.3. H II Region Star Formation Properties 

The seven radio H II regions are among the brightest 
optically in NGC 4258, with L Hq £ 10 38 erg s" 1 (Courtes 



et al. 1993). We have calculated the star formation rate, 
SFR, for massive stars from their thermal radio spectral 
luminosity, Lt, using Eq. 23 of Condon (1992): 

L T ~ 5.5 x 10 2( V- al )(SFi?) (1) 

where Lt is in W Hz -1 , v is in GHz, and SFR is in Mq 
yr _1 . Condon (1992) derives this relation in the context 
of the global star formation properties for normal galaxies 
using the "extended" Miller & Scalo (1979) initial mass 
function ip(M) oc M~ 5 / 2 , for stellar masses above 5 Mq 
and below ~100 Mq. The relatively dust-free Lt was 
determined at 6 cm assuming a distance of 7.3 Mpc to 
NGC 4258 (Herrnstein et al. 1997). For H II region candi- 
dates with —0.25 < a < —0.1 (sources 3, 4, 5, and 9), the 
uncertainties in a are significant and, therefore, we can- 
not rule out a nonthermal component. A correction was 
first applied for residual nonthermal contamination which 
was, presumably, not fully removed in the background- 
subtraction procedure discussed in Section 3. A nonther- 
mal spectral index of a nt = —0.8 was assumed for the 
contamination. This correction is at the ~1<7 level, except 
for source 9, for which the flux density is reduced by 2a 
and 3.5a at 6 and 20 cm, respectively. We present in Table 
3 for each H II region the corrected thermal flux densities 
(Sg h 20 ), Lt, the SFR, the production rate of Lyman con- 
tinuum photons (-/Vtjv, in s_1 ) ; an d the radio supernova 
rate (^sn, in yr -1 ). Ayv and z^sn were calculated using 
Eqs. 24 and 20 from Condon (1992): 

^sn ~ 0.041(5^) (2) 

A uv ~ 3.5 x 10 53 (SFR), (3) 

where the latter equation assumes that dust absorption of 
Lyman continuum photons is negligible. 

4.1.4. H II Region Luminosity Functions 

H II region luminosity functions, i.e., the number distri- 
bution of source luminosities, can reflect the mass distribu- 
tion of the parent molecular clouds and the efficiencies for 
conversion of molecular gas during star formation. Con- 
siderable foreground extinction impedes Ha studies of H II 
regions in the Galaxy. Extragalactic luminosity functions 
are important, since the lack of known distances, sample 
incompleteness, and selection biases in the Galactic sam- 
ple prevent accurate luminosity functions from being cal- 
culated. Previous radio observations of extragalactic H II 
regions, while not affected by extinction and relative dis- 
tance ambiguities, usually suffer from limited resolution 
and sensitivity and the related difficulties in nonthermal 
background subtraction. Thus far, radio luminosity func- 
tions of H II regions are available for non-Local Group 
galaxies for only M81 (Kaufman et al. 1987), M51 (van 
der Hulst et al. 1988), NGC 6946 (Hyman et al. 2000), and 
certain starburst galaxies. 

Since we have detected only seven bright radio H II re- 
gions in NGC 4258, we can only plot the brightest end of 
the luminosity function. The cumulative radio luminosity 
functions, i.e., N(S) (where N(S) is the number of sources 
with radio flux densities > S) versus S, for the seven H II 
region candidates in Table 1 are shown in Figure 4 for 
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both 6 cm (diamonds) and 20 cm (circles) wavelengths. 
The luminosity functions have nearly identical shapes and 
power-law indices (3 {N(S) oc S' 3 ), with j3 = -1.3 ± 0.2 at 
both wavelengths. 

The values we obtain for (3 are similar to those found 
in earlier studies. Recently, we reported values of f3 = 
-1.2 ± 0.2 at 6 cm and (3 = -1.3 ± 0.2 at 20 cm for the 
bright end of the luminosity function formed from 43 ra- 
dio H II region candidates detected in NGC 6946 (Hyman 
et al. 2000; Lacey et al. 1997). Viallefond & Goss (1986) 
found (3 = —1.18 ± 0.14 for the H II region 20 cm lumi- 
nosity function of M33. Kaufman et al. (1987) and van 
der Hulst et al. (1988) have fit their 20 cm luminosity 
functions for H II regions in M81 and M51 and obtained 
(3 = —1.8 ± 0.4 and (3 ~ —2, respectively, for the brightest 
sources. Israel (1980) obtained (3 values ranging from —1.1 
to —5.5 for the high end of the radio luminosity functions 
for 14 late-type galaxies, with a steeper value of —1.7 for 
NGC 4258. The difference in values probably arises from 
resolution differences and significant background contam- 
ination. We conclude that the index determined from our 
high end luminosity function is likely more reliable. We 
can clearly distinguish the H II regions in the northern 
arm from each other and can remove background contam- 
ination more accurately. 

Courtes et al. (1993) have derived the Ha luminosity 
function from measurements of 137 H II regions in NGC 
4258. However, they could not fit a single power law to the 
bright end because of a break in the slope seen at 10 38 7 
ergs s _1 (this break in a number of galaxies may be due to 
an unknown mechanism that inhibits the formation of su- 
pergiant H II regions [Kcnnicutt, Edgar, & Hodge 1989)], 
to a transition from ionization- to density-bounded regions 
[Rozas, Beckman, & Knapen 1996], or to a variable, rather 
than constant, star formation evolution [Feinstein 1997]). 
From Fig. 4 of Courtes et al. (1993) for L £ 10 38 erg s"\ 
we estimate that the slope (3 ~ — 1, somewhat shallower 
than the radio-determined slope. 

4.1.5. H II Region Extinction 

From our radio and optical measurements we have esti- 
mated the extinction at Ha, A a , toward the H II regions 
listed in Table 1, using Eq. 2 of Kaufman et al. (1987): 

A a = 2.5lo Sw [C(u,T e )S u /S Ha \ (4) 

where the electron temperature T e is assumed to be 10 
K, and C{v,T e ) = 8.01 x 10" 13 erg cm" 2 s" 1 mJy^ 1 at 
i/=1.465 GHz and C{v,T e ) = 9.05 x 10~ 13 erg cur 2 s" 1 
mJy" 1 at z/=4.885 GHz (Gebel 1968). 

We obtain mean values of 2.4± 1.0 at 6 cm and 2.5 ± 1.0 
at 20 cm for the visual extinctions, Ay, after applying a 
standard galactic reddening law, Ay — 1.28A a (Miller & 
Mathews 1972). Due to NGC 4258's high galactic lati- 
tude, the Galactic foreground extinction is negligible (Ab 
= 0.00; NED). We find no correlation between Ay and 
radio flux density, which implies that the optical and ra- 
dio H II region luminosity functions should have similar 
shapes. Such a correlation was also not found for H II 
regions in NGC 6946 (Hyman et al. 2000), M81 (Kaufman 
et al. 1987), and M51 (van der Hulst et al. 1988). 

The highest extinction is found for the combination of 
sources 5 and 6, which are located close to the center of 



NGC 4258, where more dust might be expected. We also 
have obtained an estimate of ^5.5 mag for the extinction 
toward the flat spectrum nucleus which is comparable to 
that for sources 5 and 6, and significantly greater than for 
those with lower extinctions located further from the cen- 
ter. Recently, Chary et al. (2000) from high-resolution in- 
frared observations estimate a visual extinction of 18 mag 
for the toward the compact (size < 7 pc) source at the 
center of NGC 4258. 

The model by Sarazin (1976) predicts that internal dust 
in H II regions would give rise to an extinction gradi- 
ent with galactocentric radius, related to metal abundance 
gradients. While a shallow metal abundance gradient (Du- 
til & Roy 1999) and a moderate extinction gradient (Dutil 
1998) are found from optical data alone for NGC 4258, our 
sample of extinction measurements is too limited to show 
a clear trend. However, an extinction gradient was also 
not found for NGC 6946 (Belley & Roy 1992; Hyman et 
al. 2000), M81 (Kaufman et al. 1987), and M51 (van der 
Hulst et al. 1988). 

Our mean Ay values are about 1 mag higher than the 
mean value of 1.3 mag obtained for 122 H II regions in 
NGC 4258 using Balmer decrements (Dutil & Roy 1999; 
Dutil 1998). Their extinction values, ^ly(Bal), for sources 
corresponding to our H II region candidates are listed in 
Table 1. The largest differences (from ^0.5 to ^2 mag) 
are for sources 1,9, and the combination of 5 and 6, and 
the smallest differences (<^0.2 mag) are for sources 2, 3, 
and 4. In other galaxies, such as NGC 6946 (Hyman 
et al. 2000), M33 (Viallefond & Goss 1986), M51 (van 
der Hulst et al. 1988), and M81 (Kaufman et al. 1987), 
the radio-determined extinction is also higher than the 
Balmer decrement extinction. This discrepancy has been 
attributed to patchy dust internal and external to H II 
regions, to aperture differences between the radio and op- 
tical measurements, or to the variation of line ratios over 
the spatial extent of the source. 

A comparison of extinction estimates using the radio- 
to-optical method for other galaxies shows that our mean 
value of 2.5 ± 1.0 for NGC 4258 is similar to that of 
2.2 ±0.7 for M51 (van der Hulst et al. 1988), but is greater 
than 1.0 ± 1.2 for NGC 6946 (Hyman et al. 2000) and 
1.0 ± 0.4 for M81 (Kaufman et al. 1987). 

4.2. Supernova Remnants 

The SNR candidates, sources 7, 8, and 10-12, for a dis- 
tance of 7.3 Mpc, have spectral luminosities at 6 cm of 
1.6, 1.4, 1.2, 1.7, and 2.0 x 10 22 erg s" 1 Hz" 1 , respectively. 
These luminosities are all ~2-3 times the 6 cm luminosity 
of Cas A, and are comparable to the extraordinary and 
very luminous SNR in NGC 6946 (Van Dyk et al. 1994). 
Samples of extragalactic SNRs, not subject to the extinc- 
tion and source confusion problems of Galactic SNRs, al- 
though growing, are still quite small, and these very bright 
examples are particularly intringuing, since they may the 
remnants of very massive single stars (Dunne, Gruendl, & 
Chu 2000) or may be multiple interacting SNRs (Blair, Fe- 
sen, & Schlcgcl 2001). These bright sources are deserving 
of further observational attention. No fainter likely radio 
SNRs are detected in NGC 4258. 

Cecil et al. (2000) claim that what is our source 7 is a 
putative northern radio "hot spot" created by the nuclear 
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jet impacting the ambient medium. Although we clearly 
agree that the source is nonthermal, Cecil et al. (2000) 
measure a size of 14" at l'/3 resolution at 20 cm; at the 
same frequency we measure at 4" resolution a size of at 
most 5" (FWHM). Wc find the size difference hard to rec- 
oncile and speculate on how much background emission 
was included in their measurement. We therefore call into 
question its identification as a hot spot of the jet out- 
flow. Similarly, Cecil et al. (2000) isolate a corresponding 
southern hot spot; this source is even more embedded in 
the background, and Figure la shows it to be complex 
in morphology. Therefore, we have not included it in our 
detected source list. This source is most certainly nonther- 
mal as well, with a between —0.7 and —0.4, but, again, 
the diffuse background is too bright to make an accurate 
estimate. 

4.3. Nucleus 

In our concatenated dataset, the nucleus of NGC 4258 
has a relatively flat spectrum, a = — 0.24±0.10, between 6 
and 20 cm. Turner & Ho (1994) found from observations in 
1982 and 1983 an inverted spectrum, a = +0.4, between 2 
and 6 cm, characteristic of an opaque synchrotron source. 
We undertook mapping the individual observations at each 
frequency made only in the A-configuration (to image the 
nucleus at the highest resolution and to minimize back- 
ground contamination), to look for possible variability in 
the radio emission from the nucleus. We measured the in- 
tegrated flux density of the source using the tasks IMFIT 
and IMEAN. (Uncertainties in the measurements are de- 
rived in a similar fashion to those for SN 1981K by Van 
Dyk et al. [1992].) The sampling in time is not very fre- 
quent, but Figure 5 indicates that during the monitoring 
of SN 1981K, the flux density varied by as much as 47% at 
both 6 and 20 cm, possibly undergoing outburst activity. 
The spectral index varied from —0.18 to —0.41 during this 
interval. To our knowledge, this is the first representation 
of the variability of the nucleus of NGC 4258. Variabil- 
ity of this kind is also exhibited by the LINER nucleus of 
M81 (Ho et al. 1999), albeit in much greater detail. Such 
behavior is consistent with the picture of a supermassive 
black hole being responsible for the activity in the radio 
from the nuclei of both galaxies. 

4.4. The Anomalous Radio Arms 
AAA. Description 

Our high resolution images of NGC 4258 at 6 and 20 cm 
(see Figures la and b) show greater detail in the anoma- 
lous radio arms not visible in the previous lower resolution, 
lower sensitivity maps presented elsewhere (e.g., Hummel 
et al. 1989; van Albada & van der Hulst 1982), in par- 
ticular, a large number of knots and filaments along both 
arms. As noted by previous studies, the general appear- 
ance of the northwestern and southeastern arms is quite 
different. Many of the individual sources listed by Turner 
& Ho (1994) along the arms are actually just the peak 
emission in the arms of the various knots. 

The brightness of the northwestern anomalous arm de- 
creases first to the north of the central region of NGC 
4258 and then brightens again by ^2.5 times at both 
wavelengths, to a maximum of 0.54 mJy beam -1 and 1.2 



mjy beam -1 at 6 and 20 cm, respectively. The half- 
power width of this northwestern arm at the position of 
maximum flux density before bifurcation is 7'.'2 ± 1" af- 
ter correcting for the beamsize broadening. The bright- 
ening and the subsequent bifurcation of the northwestern 
anomalous arm about ~1' from the nucleus is apparent on 
even very low resolution maps. The northern and south- 
ern subsections of this arm after bifurcation also dim and 
then brighten considerably after bifurcation. While the 
northern section widens and further fragments into several 
bright filamentary structures, the southern section remains 
narrow and brightens into only one main filament. 

In contrast to the northwestern anomalous radio arm, 
the southeastern anomalous arm has only one relatively 
bright region along its length. This region is located ^3 
times further along the arm from the center of NGC 4258 
than the bright region in the northwestern arm before bi- 
furcation. The flux density of the intervening length of the 
southeastern arm dims significantly from the central region 
to a minimum of ^0.05 mJy beam -1 at 6 cm and ^0.15 
mJy beam -1 at 20 cm, and then increases to a maximum 
of 0.42 mJy beam -1 and 0.86 mJy beam -1 at 6 and 20 
cm, respectively. The beam-corrected, half-power width 
of the southeastern arm at the position of maximum flux 
density is 6'.'7± 1". 

Emission (at the > 3a level) is detected from the 20 cm 
radio "plateaus" (see Figure 6a) south of the northwest- 
ern arm and north of the southeastern arm. We detect 
both broad, diffuse emission and a few large structures 
(> 5tr level) in the western plateau, while very little dif- 
fuse emission is evident from the eastern plateau. The 
eastern plateau does show several structures (> 3a level) 
extending into it from the southeastern arm. Since the ra- 
dio emission is largely nonthermal (see §4.4.2 below), the 
plateaus and any substructuring are much fainter at 6 cm 
(see Figure 6b). 

While much of the plateau emission is faint in these high 
resolution maps, comparison of our lower resolution maps 
(not shown) with those of others (Hummel et al. 1989; van 
Albada & van der Hulst 1982) reveals more distinct diffuse 
emission and extensions from the arms into the plateaus 
at both wavelengths. Also, the absolute flux level of the 
plateaus is underrepresented at 6 cm, and less so at 20 
cm, due to missing very short spacings in the UV cover- 
age. While D configuration observations and zero-spacing 
fluxes were included in the imaging process, a shallow (~1- 
2a), negative "bowl" is still present in the 6 cm image, and 
a shallower one in the 20 cm image, further hindering the 
detection of fainter emission in the radio plateaus. 



4.4.2. Spectral Index Map 

Hummel et al. (1989) determined the spectral index dis- 
tribution of the radio arms at 14" resolution, and de Bruyn 
(1977) and Krause et al. (1984) at even poorer resolution. 
No large-scale spectral index variations were observed, al- 
though Hummel et al. (1989) detected a possible flattening 
of the spectrum from a = —0.7 to —0.5 along the northern 
ridge of the northwestern arm. We have generated a spec- 
tral index map of the radio arms from our MEM images to 
search for variations that may be more evident at higher 
resolution and sensitivity. 
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In order to make the spectral index map, we first ad- 
justed the absolute flux density level of the images as 
corrections for the negative "bowls" discussed above in 
§4.4.1. The background level was sampled along the rel- 
atively emission-free regions below the southeastern arm 
and above the northwestern arm and determined to be —20 
± lO^Jy beam" 1 at 6 cm and —15 ± 10 /iJy beam -1 at 
20 cm. These mean values were subtracted from the cor- 
responding images. Then, the AIPS task COMB was used 
to create the spectral index map from the bowl-corrected 
images. We applied a clipping level of 90/iJy beam -1 (4.5a 
at 6 cm, 3cr at 20 cm) in order to blank out regions of low 
signal-to-noise ratio which tend to give unrealistic spectral 
index values. The resulting spectral index map is shown 
in Figure 7. 

The nucleus and H II region candidates are obvious flat 
spectrum sources with a between +0.05 and —0.25 (see 
Table 1), while the radio arms have a generally steep spec- 
trum with a = —0.65 ± 0.10. No large-scale spectral in- 
dex variations are apparent along either arm, in agreement 
with Krause et al. (1984). In addition to the spectral index 
map, in Figure 8 we present mean values of a measured 
using rectangular apertures typically 8" x 8" in size at a 
number of positions along the anomalous radio arms, in 
order to further investigate possible variations not appar- 
ent in Figure 7. The uncertainties reflect both the stan- 
dard deviation in the spectral index measurements within 
each aperture and uncertainty in the bowl correction. The 
points plotted in Figure 8 are derived from the highest ra- 
dio surface brightness regions of the anomalous arms mea- 
sured from the nucleus outward. The spectral indices for 
the northwestern and southeastern arms are quite similar; 
however, the spectrum of the northern and southern sub- 
sections of the northwestern arm appears to flatten after 
bifurcation from ~ —0.7 to ~ —0.3. No such trend is seen 
for the southeastern arm. 

Note that we are not able to study the spectral index 
distribution of the radio plateaus, due to poor signal-to- 
noisc in those regions. Using ~1' resolution images at 0.6 
and 1.4 GHz, de Bruyn (1977) determined a to be -0.7, 
consistent with the arms. For the integrated radio emis- 
sion, Krause et al. (1984) and Hummel et al. (1989) report 
a spectral index of a — —0.60 ± 0.03 for frequencies be- 
low 5 GHz, in agreement with our value for the arms, and 
— 1.16 ± 0.05 at higher frequencies. 

We do not see in either Figure 7 or 8 any significant 
variation of spectral index across the anomalous arms, as 
claimed to be the case by Hummel et al. (1989). Such a 
variation is used by Cox & Downes (1996) as an argument 
that the anomalous arms are shocked gas flow due to a 
bar potential. Any variations are on a small spatial scale 
without any systematic trends. We see nothing in our 
maps to contradict the contention that the arms are asso- 
ciated with changing jet activity from the nucleus (Cecil et 
al. 2000) and are in the galactic plane. Cecil et al. (2000) 
claim that the jet has recently moved out of the plane, but 
above we have already called into question the identity of 
their putative radio jet hot spots. 



5. CONCLUSIONS 

We have mapped the radio structure of the galaxy NGC 
4258 at 6 and 20 cm with higher resolution than previously 
available to identify H II region and SNR candidates and 
to determine the structure of its peculiar anomalous radio 
arms. Seven thermal and five nonthermal radio sources are 
detected on our 3'.' 5 resolution radio maps that we iden- 
tify as H II region and SNR candidates, respectively. We 
have identified optical counterparts to all of the H II region 
candidates and to three of the SNR candidates on the Ha 
image by Dutil & Roy (1999). Two of the SNR candidates 
also coincide with X-ray sources (Vogler & Pietsch 1999). 

By comparing the radio and optical measurements of 
the H II region candidates for the first time, we obtain an 
average value of Ay = 2.5 ± 1.0 for the visual extinction 
toward the H II regions in NGC 4258. This extinction 
in NGC 4258 is similar to that found for M51 and -1 
mag higher than that found for M81 and NGC 6946. We 
obtain for NGC 4258 higher values of extinction from our 
radio-to-Ha method than those obtained using the Balmer 
decrement (Dutil 1998), which may be due to a patchy 
dust environment within or near the H II regions. 

We find the bright end of the H II region radio luminos- 
ity function to be fit well by a power law with an index 
of P = -1.3 ± 0.2 at both 6 and 20 cm. This value is 
within the range of indices derived for other nearby spiral 
galaxies. 

The SNRs are found to be very luminous, compara- 
ble to the extraordinary SNR in NGC 6946 (Van Dyk et 
al. 1994). Additionally, the LINER nucleus is found to 
be variable in the radio, consistent with the behavior of 
M81's nucleus (Ho et al. 1999). 

The anomalous radio arms are found to have no large 
scale variation in spectral index, in agreement with Krause 
et al. (1984), although a possible flattening of the spectral 
index along the northwestern anomalous sections may be 
seen. We obtain a mean value of a — —0.65 ± 0.10 for 
the arms, consistent with that found in lower resolution 
studies. 
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Table 1 



Discrete Radio Sources in NGC 4258 a 


Source 


R.A.(J2000) 


Dcc.(J2000) 


s 6 


S20 






"20 






S(Ha) b 


A v (6) 


A v (20) 


A v (Bal) c 




(h m s) 


(° ' ") 


(mjy) 


(mjy) 
















(mag) 


(mag) 


(mag) 


1 


12 18 55.2 


47 20 35.5 


1.32 ± 0.08 


1.24 ± 


15 


+0 


05 ± 





11 


2.42 ± 


0.01 


2.2 ± 0.1 


2.0 ± 0.2 


1.81, 1.48 


2 


12 18 55.8 


47 20 24.3 


0.88 ± 0.08 


1.00 ± 


15 


-0 


11 ± 





15 


2.14 ± 


0.01 


1.8 ± 0.1 


1.8 ± 0.2 


1.65 


3 


12 18 56.5 


47 20 14.4 


0.64 ± 0.11 


0.81 ± 


13 


-0 


20 ± 





20 


1.36 ± 


0.01 


2.0 ± 0.2 


2.2 ± 0.2 


1.98 


4 


12 18 57.4 


47 20 04.6 


0.50 ± 0.08 


0.66 ± 


11 


-0 


22 ± 





19 


1.69 ± 


0.01 


1.4 ± 0.2 


1.6 ± 0.2 


1.33, 1.33 


5 


12 18 58.2 


47 19 03.9 


0.56 ± 0.12 


0.66 ± 


18 


-0 


13 ± 





28 












6 


12 18 58.2 


47 19 10.1 


0.33 ± 0.09 


0.33 ± 


18 


+0 


02 ± 





51 












5+6 d 






1.18 ± 0.12 


1.65 ± 


28 


-0 


26 ± 





13 


0.52 ± 


0.01 


4.2 ± 0.1 


4.5 ± 0.2 


2.26 


7 


12 18 57.0 


47 19 03.9 


0.26 ± 0.08 


0.71 ± 


14 


-0 


82 ± 





30 












8 


12 18 47.9 


47 18 03.3 


0.23 ± 0.06 


0.43 ± 


12 


-0 


50 ± 





32 












9 


12 19 01.3 


47 15 25.0 


1.89 ± 0.15 


2.56 ± 


22 


-0 


25 ± 





10 


2.58 ± 


0.02 


2.6 ± 0.1 


2.9 ± 0.1 


1.42 


10 


12 18 55.3 


47 16 47.9 


0.20 ± 0.10 


0.31 ± 


07 


-0 


37 ± 





46 


0.29 ± 


0.01 








11 


12 18 57.6 


47 16 07.0 


0.27 ± 0.09 


0.53 ± 


08 


-0 


56 ± 





31 


0.15 ± 


0.01 








12 


12 18 56.3 


47 16 50.3 


0.31 ± 0.09 


0.64 ± 


12 


-0 


58 ± 





29 


0.55 ± 


0.01 








nucleus 


12 18 57.5 


47 18 14.5 


2.36 ± 0.16 


3.14 ± 


29 


-0 


24 ± 





10 


0.46 ± 


0.01 









a The optical data arc derived from the Hct image by Dutil & Roy (1999). 
b In units of 10 -13 ergs cm -2 s _1 . 

c Extinctions determined from Balmer decrements are from Dutil (1998). For apparently single radio sources corresponding to two optical 
sources, two extinction values arc listed. 

d The Ha flux of a single optical source overlapping sources 5 and 6 was measured. See text in §4.1.1 



Table 2 



Background Radio Sources 


Source 


R.A 


(J2000) 


Dec.(,T2000) 


Se 


•520 


"20 








(h m s) 


(° ' ") 


(mJy) 


(mJy) 








13 


12 


18 56.5 


47 23 53.4 


0.67 ± 0.10 


1.81 ± 0.09 


-0.82 ± 





13 


14 


12 


18 46.1 


47 23 12.4 


0.36 ± 0.10 


0.83 ± 0.12 


-0.70 ± 





25 


15 


12 


19 19.0 


47 19 07.6 


0.92 ± 0.12 


2.09 ± 0.23 


-0.69 ± 





14 


16a 


12 


19 32.6 


47 18 20.5 


2.43 ± 0.20 


7.76 ± 0.40 


-0.97 ± 





08 


16b 


12 


19 32.0 


47 17 55.8 


1.75 ± 0.14 


2.09 ± 0.23 


-0.15 ± 





11 


16c 


12 


19 31.5 


47 17 48.4 


0.66 ± 0.20 


2.48 ± 0.11 


-1.09 ± 





26 


16d 


12 


19 31.0 


47 17 38.5 


a 


3.16 ± 0.16 








17a 


12 


18 46.6 


47 14 44.0 


1.42 ± 0.15 


5.17 ± 0.35 


-1.07 ± 





10 


17b 


12 


18 46.4 


47 14 27.9 


0.67 ± 0.17 


b 








17c 


12 


18 46.2 


47 14 19.3 


3.06 ± 0.22 


b 








17b+c 








3.73 ± 0.28 


11.68 ± 0.34 


-0.95 ± 





07 



a Source 16d is not detected at 6 cm, probably due to primary beam attenuation. 

Source 17b cannot be separated from 17c at 20 cm, due to bandwidth smearing. The combined flux at each wavelength was used to calculate an 
average spectral index given below for source "17b-f-c". 



Table 3 
Star Formation Properties 3 - 



Source 


c th 


Crth 

D 20 


Z/T 


SFR 


Nuv 


fSN 




(mJy) 


(mjy) 


(W Hz" 1 x 10 1S ) 


(M Q yr" 1 ) 


(s- 1 x 10 51 ) 


(yr" 1 X 10- 4 ) 


1 


1.32 


1.24 


8.44 


0.018 


6.29 


7.37 


2 


0.88 


1.00 


5.61 


0.012 


4.18 


4.90 


3 


0.58 


0.65 


3.70 


0.008 


2.76 


3.23 


4 


0.45 


0.50 


2.84 


0.006 


2.12 


2.45 


5 


0.55 


0.62 


3.49 


0.007 


2.61 


3.05 


6 


0.33 


0.33 


2.14 


0.005 


1.59 


1.87 


9 


1.60 


1.80 


10.2 


0.022 


7.61 


8.91 



a Star formation properties were calculated based on 5g h and equations in Condon (1992). For sources with a < —0.1, the flux densities given in 
Table 1 (S6 and S20) arc corrected for residual nonthermal contamination (§4.1.3) to obtain the thermal components (5g h and S^q) listed here. For 
sources with a > —0.1, 5g h and S^q are set equal to the values of S6 and S20- 
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Fig. 1. — (a) 20 cm 3'.'5 resolution image of NGC 4258 made using the Maximum Entropy Method. The labels refer to 
sources discussed in §4.1.1 and 4.1.2. 



Fig. 1. — (6) 6 cm 3'.'5 resolution image of NGC 4258 made using the Maximum Entropy Method. 



Fig. 2. — Contour map of the Ha image from Dutil & Roy (1999) overlaid on our 20 cm 3'.'5 resolution greyscale image 
of NGC 4258 (from Figure la). The Ha contour levels are 7.5xl0~ 18 x (3, 5, 10, 25, 50, 175, 300, 500) ergs cm" 2 s" 1 . 
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12 19 02 00 18 58 56 54 52 

RIGHT ASCENSION (J2000) 

Fig. 3. — Contour map of the southern arm region of the Ha image from Dutil & Roy (1999) overlaid on the corresponding 
region of our 20 cm 3'.'5 resolution greyscale image of NGC 4258 (from Figure la). Note that all pixels with flux densities 
less than 90 mJy beam -1 (3<r) have been set to zero in order to better depict the radio source emission in the southern 
arm. The numbers on the figure correspond to sources discussed in §4.1.1. The centers and sizes of the boxes near sources 
11 and 12 indicate, respectively, the positions and positional uncertainties of two X-ray sources detected by Voglcr & 
Pietsch (1999). The Ha contour levels are 7.5xl0 -18 x (20, 50, 100) ergs cm -2 s" 1 . 
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Fig. 4. — Cumulative H II region luminosity functions at 20 cm (circles) and 6 cm (diamonds) for NGC 4258. Shown are 
fits to the bright ends of the luminosity functions at 20 cm (solid line) and 6 cm (dashed line), and have power-law index 
(3 = — 1.3 ± 0.2 at both wavelengths. 
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Fig. 5. — Flux density at 6 cm (open squares) and 20 cm (filled triangles) of the nucleus of NGC 4258, derived from 
A-configuration (highest resolution) VLA observations, as a function of time (relative to the optical maximum of SN 
1981K in the galaxy; see Van Dyk et al. 1992). The plot is suggestive of long-term variability from the nuclear source. 

Fig. 6. — (a) 20 cm 3'.'5 resolution contour map of NGC 4258, emphasizing the emission from the radio plateaus. The 
contour levels are 30 /zJy beam" 1 x (3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 35, 40, 45, 50). 



Fig. 6. — (b) 6 cm 3'/5 resolution contour map of NGC 4258. Due to its nonthermal spectrum, less of the broad disk 
emission is detected. The contour levels arc 20 ^Jy beam" 1 x (3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 35, 40, 45, 
50). 



Fig. 7. — Map of the spectral index a (S cx v +a ) between 6 and 20 cm for NGC 4258. The range of spectral index values 
are approximately indicated by the greyscale bar shown at the top of the figure. 
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Fig. 8. — Spectral index distribution along the anomalous radio arms. The circles are samples of the southeastern arm 
(negative distances). The squares represent the northwestern arm before bifurcation and the southern subsection of the 
northwestern arm after bifurcation; the diamonds represent samples of the northern subsection of the northwestern arm 
after bifurcation (positive distances). The distribution is essentially uniform, although the northwestern arm may flatten 
slightly after bifurcation. The flat spectrum point at a = —0.22 corresponds to the nucleus. 
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